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We show that in the Landau gauge of the SU(2) Yang-Mihs theory the resid- 
ual global symmetry supports existence of the topological vortices which resemble 
disclination defects in the nematic liquid crystals and the Alice (half-quantum) vor- 
• tices in the superfluid ^He in the A-phase. The theory also possesses half-integer 

. and integer-charged monopoles which are analogous to the point-like defects in the 

I nematic crystal and in the liquid helium. We argue that the deconfinement phase 
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transition in the Yang-Mills theory in the Landau gauge is associated with the pro- 
liferation of these vortices and/or monopoles. The disorder caused by these defects 
is suggested to be responsible for the confinement of quarks in the low-temperature 



^ . phase. 

PACS numbers: 12.38.-t,12.38.Aw,61.30.-v 

^ ■ The confinement of color in Quantum Chromodynamics is one of the unsolved problems 

! in the quantum field theory. Nowadays the commonly accepted point of view is that con- 
O I finement is encoded in the dynamics of gluon fields which are described by the Yang-Mills 
(YM) theory. Popular approaches to the color confinement are based on the dual supercon- 
ductor mechanism and on the center vortex picture These approaches suggest that 
the confinement is caused by a particular dynamics of special gluon configurations called 
the Abelian monopoles and the central vortices, respectively. Besides the different dimen- 
sionality of these configurations, the formulation of these confinement mechanisms requires 
a gauge fixing of the non-Abelian symmetry which is also different for both approaches. 
Recently, an interest to describe the confinement mechanism in the Landau gauge in terms 
X ; of the particular configurations of the gluon fields has emerged 0, 0. The Landau gauge 
' is very attractive because it is very well studied both in perturbative and non-perturbative 
approaches. Moreover, this gauge is color-symmetric and it can be well formulated (up to 
the Gribov-copy problem with most of other gauges also have) both in the perturbative 
regime and in the continuum space-time. 

In Ref. the dual Meissner effect caused by magnetic displacement currents was observed 
numerically in the Landau gauge. The dual Meissner effect leads to the confinement of quarks 
via formation of the chromoelectric string. This string is similar the Abrikosov vortex [5(| in 
the ordinary superconductors. 

Another approach utilizes the spin-charge separation idea 0, 0] which is widely used 
in condensed matter physics to describe the behavior of strongly correlated electrons in 
high-Tc cuprate superconductors In Ref. 0] the YM theory in the Landau gauge was 
re-formulated as a nematic liquid crystal in an internal space-time. In this paper we observe 
a relation between the YM theory and a nematic crystal without the use of spin-charge 
separation of the gluon field. 

We show that the existence of a residual global symmetry - which remains intact after the 
conventional Landau gauge fixing - leads to existence of specific topological defects which 
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are analogous to the defects in the nematic liquid crystals and in the superfluid ^He. There 
are many examples of physical systems in which the global symmetries play an essential 
role: ferromagnetics, liquid Hca and He4, various liquid crystals etc. We point out that the 
(residual) global gauge symmetry in the Landau gauge is a symmetry which can not be 
simply ignored. 

We are working with the pure SU(2) Yang-Mills theory in the Euclidean space-time. The 
Landau gauge is defined as a minimization of the gauge fixing functional, 

mmF[A^], F[A]= j d'x[Al{x)\\ (1) 

over the gauge transformations Vt G 5'0(3)gaugo- The corresponding local differential gauge 
condition is d^A^^ = 0. The condition ([T)) fixes the S'O (3)gauge gcLuge color freedom up to the 
S'0(3)giobai global color group, 

•S'O (3) gauge 5'0 (3) global , (2) 

because the gauge-fixing functional (|H) is invariant under the global (coordinate-indepen- 
dent) transformation, 

Alix) fi^f Alix) , figi G 50(3)giobai . (3) 

Since the gluon field does not transform under the action of the center Z12 of the gauge 
group SU{2), the gauge symmetry of the pure SU(2) Yang-Mills theory is in fact 5*0(3) ~ 

Sf/(2)/^2- 

The vector gauge field transforms in the fundamental representation of the residual 
S'O (3)giobai group, Eq. Q. One can construct the composite "color-spin" field 

C^\x) = A^(a;)A^(x) ■ 6^^' - A';{x)A'^{x) , (4) 

which is a scalar with respect to space-time rotations and a rank-2 tensor with respect to 
the global color transformations^. The field C transforms in the adjoint representation of 
the residual S'O (3) global group, 

C{x) ^QC{x)Q^ , (5) 

where the superscript T stands for the transposition in the color indices. By construction, 
the field C is symmetric, C = C^. 

At this point it is convenient to introduce a mechanical interpretation of the matrix O"*. 
The definition (0]) is similar to the moment of inertia tensor of a solid body composed of 
mass-centers (labeled by the integer n) with the masses located at the positions r„: 

r^ = ^m„(r^5^^-rX). (6) 

n 

At each point x of the space-time the matrix O"'', Eq. (j^, corresponds to the inertia tensor 
of a "solid body" consisting of four (/i = 1, . . . , 4 in the four-dimensional space-time) mass- 
centers with equal "masses" mi = . . . = 1. The gauge field plays a role of a coordinate 



Note that a first intension to construct the pure octet representations in the form A'^A'^^ — (5"^ 
instead of Eq. |0J would make a mechanical analogy discussed below a bit obscure. 
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of the /i**^ mass-center. The color transformation (0) of the matrix C corresponds to a 
transformation of the inertia tensor © under a spatial rotation. 

At each point x of the space-time the matrix C°'^{x) can be diagonalized with the help 
of the 5*0(3) transformation 6(x), 



The eigenvalues c^, k = 1,2,3 can be interpreted as "moments of inertia" defined with 
respect to the orthonormal "principal axes of inertia" e^, k = 1,2, 3, respectively. The axes 
Bk are normalized eigenvectors of the "moment of inertia tensor" C. The transformation 
Q{x) relates the default basis in the color space with the basis of the principal axes of inertia 
at the point x. The eigenvalues of the matrix C are ordered, ci ^ C2 ^ C3, and below we 
assume that these eigenvalues are not degenerate unless it is stated otherwise. 

To pursue this mechanical analogy further, the tensor C can be associated an "ellipsoid 
of inertia" in the color space. The directions of the axes of this ellipsoid coincide with 
principal axes of inertia e^, k = 1,2, 3, while the axes lengths are proportional to the inverse 
square roots of the corresponding moments of inertia. The direction of the longest axis of 
the ellipsoid of inertia - corresponding to the lowest principal moment of inertia - defines a 
unit vector n(x) = 63 (x) in the color space. 

The vector n(x) is direction-less since the vectors ±n correspond to the same ellipsoid of 
inertia. Indeed, the ellipsoid is invariant under the transformations of the D2 group which 
is a finite subgroup of the SO (3) group of rotations. The elements of D2 are vr-rotations 
about any two principal axis of the ellipsoid. The product of these elements generate the 
vr-rotation about the third axis. The vr-rotation of the ellipsoid about either ei or 62 axis 
leads to the flip n — n. Thus, the eigenvector corresponding to the lowest eigenvalue of 
the composite symmetric field Q defines an arrowless vector in the color space n - or, in 
other words, a projective plane RP(2) element - for a given configuration of the gauge fields 

Aa 

The arrowless feature of the vector n(x) allows to formulate a liquid crystal-like structure 
in the YM fields. The spin field n may also be associated with the largest principal axes of 
the axially symmetric molecules {i.e., which is also called the "director field") in nematic 
liquid crystals ^. The ordinary nematic crystals 0| are liquids composed of rod-like 
molecules which are randomly oriented in the isotropic (high-temperature) phase. At low 
temperatures the systems goes into the liquid crystal phase in which the "rods" tend to 
align parallel to the direction no in macroscopic volumes. The molecules in liquid crystals 
do not have a positional order while being orientationally ordered. 

The rod-like molecule in the nematic is invariant under (i) the ^2 group consisting of 
the TT-rotations about any axis perpendicular to the longest axis of the molecule and (ii) 
the 5*0(2) group of rotations about the longest axis. Thus, the physical space of the axial 
molecule - corresponding to the ellipsoid of inertia defined by the color tensor Q - is the 
coset 



C''\x) = e{x) diag(ci,C2,C3)e^(x). 



(7) 



G/H = SO{3)/{^2 X S0{2)) . 



(8) 



^ Note that in the hquid crystals the global rotations of the director are related to the rotations of the 
coordinate space, while in the YM model they are independent. However, this difference is inessential and 
it does not influence the topological and dynamical statements made below. 
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Let us imagine for a moment that we have integrated out in the YM partition function (fixed 
to the Landau gauge) all the degrees of freedom but the spins n. Then at each point of 
the space-time the physical group of the n-spins would become exactly the coset group (jHI). 
Thus, after all degrees of freedom but n are integrated out, the D2 invariance group of the 
"ellipsoid of inertia" must be replaced by the ^2 x S0{2) group which leaves the arrowless 
vector n intact. 

One should note that the coset group (jH)) does not correspond to the symmetry breaking 
G ^ H m the YM theory. In the pure YM theory the direction of the vector field n changes 
from point to point leading to the vanishing vacuum expectation value, (n) = since the 
color symmetry is obviously not broken in this case. Equation (jH)) defines the group, the 
non-unit elements of which make the physically distinguishable changes to the director field 
n. 

Despite the absence of the symmetry breaking, still one can define defects associated with 
topological invariants of the arrowless vector n. These definitions do not depends on the 
dynamics of the model, while they do depend on the possibility to define a unit vector with 
particular symmetry properties (jHI). The coincidence of the physical group (jHI) with the one 



of the nematic crystals |9|, |ll| allows us to discuss topological (with respect to the residual 
50 (3) global group) defects in the Landau gauge of the YM theory using the corresponding 
condensed matter analogues. 

The physically interesting are the first four homotopy groups of the physical space (jHj). 
The ttq group is trivial, 71q{G/H) = 1, teUing us that there are no topologically stable domain 
walls made of the director field n. 

The first homotopy group is nontrivial, 'Ki{G/H) = indicating that there are topo- 
logically stable vortices characterized be the only nontrivial element of the group. The 
^2- vortices are very well known as "disclination defects" in the physics of liquid crystals j^. 
The typical example of the straight static vortex going perpendicularly to the x — y plane 
is described by the director field of the form: 

Hvort = (cos rup, sin rup, 0) (9) 

where is the azimuthal angle in the x — y plane and n is the topological number which 
may take integer and half-integer values. The integer-values vortices are unstable, and 
the topologically stable elementary vortex has the fractional circulation number n = 1/2. 
The vortices with n = 1/2, n = 1 and n = 3/2 are visualized in Figures ^a),(b) and 
(c), respectively. In the center of the vortex the color tensor (@)) must be at last double 
degenerate to support the two-dimensional hedgehog-like singularity (jHI). 

The presence of the half-winding vortices is possible because the directions n and — n are 
physically indistinguishable. As it can be seen from Figure ^a) the director field n rotates 
in 180° degrees in the color space, {i.e., n — n), as one makes the full 360°-turn around 
the vortex position. These vortices share similarities with the central vortices observed 
numerically in the Maximal Center gauge of the YM theory, as well as with the so-called 
Alice strings discussed both in the particle physics 2| and in the condensed matter physics 
{i.e., as half-quantum vortices which exist both in nematic liquid crystals E3| and in 
A-phase of the superfluid '^He 1^ l3|). 

The non-triviality of the second homotopy group, tt2{G/H) = with = Z/Z2 = 
0,1,2,... guarantees the presence of the topologically stable monopoles. The monopoles 
are particle-like objects characterized by positive half-integer winding numbers m = M/2, 
M G which count the number of sweeping of the spatial sphere into the sphere of the 
color group. Due to the ^2^dentification n — > — n, the number m can be half-integer. 
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FIG. 1: The schematic behavior of the vortex director field nvort; Eq. ©, in the color space. 
Vortices with (a) n = 1/2 (topologically stable), (b) n = 1 (unstable) and (c) n = 3/2 (unstable) 
winding numbers are shown. 



The monopole current can be formulated in terms of the winding number, 

1 dn°- dn^ dn'^ 

^ " 8^^^^^^^"'"'^ dx dx dxR ■ 

The current is represented as 5-singularities on a set of closed loops C parameterized by the 
four- vector X(r), 



kf,{x) 



J2mc fdr?^5^'\x-X'ir)), (11) 



where mc is the charge ascribed to the loop C. By construction ()10|1 . the monopole current 
is closed, d^k^ = 0. A three-dimensional analogue of Eq. (fTIH) was used in Ref. [l^ to 
formulate 7r2-singularities in the superfluid Helium and in the liquid crystals. Note that at 
the center of the monopole the color tensor Q must be triple-degenerate to support the 
three-dimensional hedgehog-like singularity (fT^ . 

Due to the identification n ^ — n the monopole winding numbers ±mc correspond to 
the same monopole according to Eq. (fTUI) . This feature is perfectly consistent with the 
presence of the Alice .^2-vortices discussed above, because the monopole with the charge 
m transforms into the (physically the same!) monopole with the charge — m after circling 
around the Alice string |jl4| . 

A class of the static unit director fields corresponding to the m-charged monopoles is 
given m Refs. [il[l3: 



rimoii = (sin ■i?(^, m) ■ cos rmp, sin ^9(6', m) ■ sin rmp, cos rmp) , (12) 

where 9 and ip are polar and azimuthal angles, and d{9,m) = 2 arctan {[tan(^/2)]'^}. The 
positive topological number m may take integer and half-integer values. The current (jllll 
corresponding to the configuration (fT^ can be calculated with the help of Eq. (jTU)), A;^(x, t) = 
'm S . Note that the director field for the m = 1 monopole is just the standard 
hedgehog, rimon = x. 

The examples of the static m = 1/2 and m = 1 monopoles (fT^ are visualized in a 
spatial time-slice in Figures |2ta), (b), respectively. The bottom projection of the m = 1/2 
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FIG. 2: The same as in Fig. ^but for monopoles (|12p with the winding numbers (a) m = 1/2 and 
(b) m = 1. The projections of the director field n into two planes are also shown. 



monopole, Figure Efa), shows a shadow of the attached = 1/2 vortex, Figure ^^a). 

The third homotopy group is also non-trivial, Ti^iG/H) = which leads to the appear- 
ance of the instantons in the four dimensional sense. 

As an example of a particular configuration of the gluon fields, consider the BPS monopole 
solution to the SU(2) Yang-Mills equation of motion: 

At = - fir)e,atx' , /(r) = - (l - , (13) 

g r \ smh r / 

Al = -h(r)x^, /i(r) = -f^cothr- iV (14) 

g r V / 

where = a;"/|x|, and r = \x\ is assumed to be scaled by an arbitrary factor, r — > ro ■ r, 
to make a dimensionless quantity. This static configuration satisfies the local Landau gauge 
condition, = 0. 

The color tensor C"", Eq. (0]), evaluated at the monopole configuration ()13I14|) has three 
eigenvalues, Ci = C2 = {P{r) + h?{r))/g'^ and C3 = 2p{r)/g'^. Since the two out of the 
three eigenvalues are degenerate, the "ellipsoid of inertia" is axially symmetric. Moreover, 
/(r) < h{r) for r > 0, therefore Ci(r) = C2(r) > c^{r), and the ellipsoid is always rod-like 
(contrasted to a disk- like form). The longest axis of the ellipsoid forms a hedgehog-like 
structure around the center of the monopole, n = es = x. Thus, in the Landau gauge the 
BPS YM- monopole can be identified with the m = 1 charge global monopole (fT^ . which 
can be called as "arrowless hedgehog" due to the ^2 symmetry n ^ — n. Note that in the 
center of the monopole the tensor (jl)) is triple-degenerate, Ci(0) = 02(0) = 03(0) = 0. 

Before discussing the dynamics of the liquid crystal defects in the YM theory let us men- 
tion an attempt (18| to construct an effective nematic theory from the YM theory in an ex- 
plicitly gauge-invariant form. The nematic variables were constructed from the field-strength 
tensor of the SU(2) gauge field. The numerical results [3| have revealed a dominance of the 
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ultraviolet fluctuations in the dynamics of the SU(2) gauge- invariant nematic variables both 
in confinement (low temperature) and in the deconfinement (high temperature) phases. We 
believe that the gauge-fixing used in our approach smoothes the gauge fields making the 
nematic field less random. 

Since the color symmetry is not broken in the YM theory the director field does not 
have a constant direction in the color space, (n) = 0. Therefore both the confinement and 
the deconfinement phases correspond - in the language of the nematic liquid crystals - to 
the (color) isotropic phases, where the density of the topological defects is non-zero due to 
large-angle fluctuations of the director field. 

In the confinement phase a network of the defects is likely to be percolating so that any two 
infinitely-separated points have a finite probability to be connected by a world-trajectory 
of the defects. Generally, the percolation of a defect trajectory means an existence of a 
non-vanishing low-momentum (infrared) component of a corresponding field, which implies 
the condensation of the defect. The percolation (condensation) feature is shared by many 
realistic theories possessing the topological defects. The percolation happens in the phases 
where the symmetry group (responsible for the appearance these defects) is unbroken. The 
relevant examples are the symmetric phases of the Abelian Higgs (or, the Ginzburg-Landau) 
model [l^ and the non-Abelian Higgs (or, the Standard Electroweak) model where 
the monopole-like and vortex-like defects are known to be proliferating at infinitely large 
distances. Note that in nematic crystals the monopole density (contrary to the vortex 
density) is expected to be low due to a topological argument of Ref. |2^. Nevertheless, a 
low-density monopole network can also be percolating. 

To discuss the behavior of the monopoles in the deconfinement phase we note, that both 
integer and half-integer monopoles in the Landau gauge are characterized by the specific 
(hedgehog- like, as in the examples shown in Figure 0) behavior of the director field in a local 
three-dimensional time-slice perpendicular to the direction of the monopole current (fTT|) . 
In order to support a non-zero sphere-to-sphere winding number ()10|) the director field 
must evolve in all three dimensions in the local vicinity of the monopole. The hedgehog 
structures involving the Euclidean time direction are suppressed at high temperatures by 
high Matsubara frequencies, leading to a suppression of the monopoles moving in spatial 
directions. 

Thus, at high temperatures the monopoles must be static, and, as a consequence, non- 
percolating. The high and the low temperature regimes must be separated by a percolation 
transition, which is very likely to at the same temperature as the deconfinement phase 
transition as it happens in the the Abelian and non-Abelian Higgs theories [19, 20]. A 
similar property is also observed numerically for the Abelian monopoles in the Maximal 
Abelian gauge of the Yang-Mills theory [j^ . 

Similar arguments can be applied to the half-quantum vortices in the Landau gauge: 
they are expected to be proliferating in the confinement phase and tend to be static in the 
deconfinement phase. Since the vortices are line-like defects, two different possibilities can 
be in the deconfinement phase: (i) the vortices may exist in the form of dominantly static 
and non-percolating short loops, (ii) or the static vortices may exist in the form of long loops 
leading to a spatial percolation of the vortex trajectories. The experience gained with the 
center vortices in the Maximal Center gauge ^] tells us that it is the first option that is 
likely to be realized in the high-temperature phase. 

The dynamics of the topological defects is important because it is related to the con- 
finement of color. In the confinement phase the percolating defects should cause a disorder 
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which should lead to the area law of sufficiently large Wilson loops. In the deconfinement 
phase the defects are expected to be dominantly static and the disorder spreads only in 
spatial dimensions leading to the deconfinement of the static quarks and to confinement of 
the "spatial" quarks in a sense of the area-law for the spatial Wilson loops. 

Let us discuss the (dis)order parameter phase transition within the described (nematic) 
picture. In the language of the topological defects, the disorder parameter for the decon- 
finement phase transition is given by the percolation probability for the monopoles (and, 
probably, vortices) to proliferate for infinitely long distances. This probability should be 
zero in the deconfinement phase and non-zero in the confinement phase. 

In the language of fields a best candidate to discuss the transition might be the color 
tensor C"*, Eq. because the field C"* (up to an inessential color-singlet part) looks 
similar to the diamagnetic susceptibility in the nematic liquid crystals. The diamagnetic 
susceptibility is known to be the order parameter of the nematic-isotropic phase transition 
in the liquid crystals j^. In the nematic (low temperature) phase the director field is 
dominantly constant leading to an anisotropy in the diamagnetic susceptibility, while in the 
high temperature phase the director field is random and the diamagnetic susceptibility is 
isotropic. Since the color symmetry is unbroken in the YM theory, one must have {C""^) oc 5°'^ 
both in the confinement and in the deconfinement phases. Consequently, the field C""^ can 
not be used as an order parameter for the deconfinement phase transition. 

Thus, the deconfinement phase transition in the YM must be associated with a transition 
from one isotropic phase of the nematic crystal to another isotropic phase. The two distinct 
isotropic phases were indeed observed in the lattice numerical simulations of the three- 
dimensional nematic crystals [i^. One of such isotropic phases is characterized the defect 
condensation. In the language of the YM theory it corresponds to the confinement phase. 
Another isotropic phase - called the topologically ordered phase ji^ - is characterized by 
the absence of the condensate (the analog of the deconfinement phase in the YM theory). 
Surprisingly, the isotropic-to-isotropic transition in the nematic crystal is the second-order 
phase transition lying in the 3D Ising universality class [2^ . The fact that the order and the 
universality class of the phase transition in the nematic crystal coincide with the order and 
the class of the finite-temperature phase transition in the SU(2) YM theory [2^ provides an 
an additional support in favor of the liquid crystal interpretation of the YM theory in the 
Landau gauge. 
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